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ABSTRACT 

We construct a simple accretion model of a rotating gas sphere onto a 
Schwarzschild black hole. We show how to build analytic solutions in terms 
of Jacobi elliptic functions. This construction represents a gene ral relativisti c 



generalisation of the Newtonian accretion model first proposed by lUlrichI (1l976l ). 
In exactly the same form as it occurs for the Newtonian case, the flow naturally 
predicts the existence of an equatorial rotating accretion disc about the hole. 
However, the radius of the disc increases monotonically without limit as the flow 
reaches its minimum allowed angular momentum for this particular model. 

Subject headings: accretion - relativity - hydrodynamics 



Introduction 



Steady spherically symm etric a c cretio n onto a central gravitational potential (e.g. a 
star) was first investigated by iBondil (Il952l ). This pioneerin g work turned out t o have many 
different applications to astrophysical phenomena (see e.g. iFrank et al.l (120021 )). despite of 
the fact that it was only constructed for curiosity, rather than a realistic idea to a particular 
astrop hysical situati on (IBondil I20051) . A general relativistic generalisation to the work of 
Bondi was made by Michel ( 19721 ). Both models can be seen as astrophysical examples of 
transonic flows that naturally occur in the Universe. 

Realistic models of spherical accretion require an extra ingredient that seems inevitable 
in many astronomical situations. This is so because gas clouds, where compact objects are 
embedded, have a certain degree of rotation. This rotation enables the formation of an 
equatorial accretion disc for which gas particles rotate about the central object. The first 
steady accretion model, in which a rotating gas sphere with infinite extent is accreted to 
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a central object was first investigated by lUlrichI (jl976l ). In his model, lUlrichI considered a 
gas cloud rotating as a rigid body and took no account of pressure gradients associated to 
the infalling gas. In other words, his analysis is approximately ballistic. For, the the initial 
specific angular momentum of an infalling particle is small and heating by radiation as well 
as viscosity effects are negligible. In addition, pressure gradients and internal energy changes 
along the streamlines of a supersonic flow provide negligible contributions to the momentum 



and e nergy balances respectively (cf. IUlrichlll976l : ICassen fc Moosmanlll98ll : iMendoza et al. 
200l). 



A flrst order general relatiyistic approximation of a rotating gas sphere was made by 
Beloborodov fc lUarionovl (120011 ). In their model, they used approximate solutions for the 
integration of the geodesic equation and their boundary conditions are such that the speciflc 
angular momentum for a single particle h < 2rg, where Tg is the Schwarzschild radius. 
In here and in what follows we use a system of units for which G = c = 1, where G 
is the gravitational constant and c the spe ed of light. In this article we show that such 
a model is not a general relativistic lUlrichI flow, sinc e its appropriate generalis ation must 
satisfy the inequality h > 2rg. A pseudo-Newtonian iPaczynsky fc Wiital (19801) numerical 



appro ximation of the extreme hyperbolic h = 2r^ case was discussed by lLee fc Ramirez- Ruiz 
(120061 ) . We show that this pseudo-Newtonian numerical approach differs in a signiflcant way 
when compared with the complete general relativistic solution. 

In this article we develop a full general relativistic model of a rotating gas sphere of 
inflnite extent that accretes matter onto a centrally symmetric Schwarzschild space-time. 
We assume that heating by radiation and viscosity effects are small so that the flow can be 
treated as ideal. Since pressure gradients and internal energy changes along the streamlines 
of a supersonic flow provide negligible contributions to the momentum and energy balances 
respectively, the flow is well approximated by ballistic trajectories. We also assume that the 
self-gravity of the accreting gas does not change the structure of the Schwarzschild space- 
time. This is of course true if the mass of the central object that shapes the space-time is 
much greater than the mass of the rotating cloud. With these assumptions, we find velocity 
and particle number density fields as well as the streamlines of the fiow in an exact analytic 
form using Jacobi Elliptic functions. The remaining thermodynamic quantities are easily 
found by assuming a polytropic fiow, for which the pressure is proportional to a power of 
the particle number density (see for e.g. IStanyukovichI Il960l ). In section [2] we state the 
main results from general relativity used to solve the model introduced in section [31 We 
show in secti on H] that the general solution converges to the accretion model considered by 
UlrichI (119761 ) and that for the case of a null v alue for the sp ecific angular momentum, the 
velocity field converges to the one described by iMichell (Il972l ) for a null value of the pressure 
gradients on the fiuid. The particular case of a minimum specific angular momentum h = 2r„ 
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is calculated in section [5l and it is shown that the solutions can be found with the aid of simple 
hyperbolic functions. Finally, in section [6] we discuss the physical consequences implied by 
this general relativistic model. 



2. Background in celestial mechanics for general relativity 

The main results from relativistic gravity to be used throughout the arti cle are stated 
in this se ction. The re a der is r eferred to the gene r al relat iv ity textbooks by iMisner et al. 



(|l973h: IChandrasekhad (119831 ): iLandau fc Lifshitzl (!1994ah : iNovikov fc Frolovl (119901 ) and 
Waldl (119841 ) for further details. 



It is well known that the vacuum Schwarzschild solution describing the final product 
of gravitational collapse contains a singularity which is hidden by a horizon. The solution 
corresponding to an exterior gravitational field of static, spherically symmetric body is given 
by the Schwarzschild metric: 



2M 



dt^ + 1 - 



2M 



-1 



dr^ + rMO^ 



where dfl^ = d9^ + sin^6'dv9^ represents the square of an angular displacement. The total 
mass of the Schwarzschild field is represented by M. The temporal, radial, polar and az- 
imuthal coordinates are represented respectively by t, r, 6 and (p. In equation (Q, we have 
chosen a signature (— , +, +, +) for the metric. In what follows, Greek indices such as a, f3, 
etc., are used to denote space-time components, taking values 0, 1, 2 and 3. 



Birkhoffl (119231 ) showed that it is possible to solve the vacuum Einstein field equations for 



a general spherically symmetric space-time, without the static field assumption. It follows 
from his calculations that the Schwarzschild solution remains the only solution of this more 
general space-time. 

The behaviour of light rays and test bodies in the exterior gravitational field of a spher- 
ical body is described by analysing both, timelike and null geodesies. In order to do that, 
we first note that the Schwarzschild metric has a parity reflection symmetry, i.e. the trans- 
formation 6 IT — 6 leaves the metric unchanged. Under these considerations it follows that 
if the initial position and tangent vector of a geodesic lies in the equatorial plane 6 = tt/2, 
then the entire geodesic must lie in that particular plane. Every geodesic can be brought to 
an initially equatorial plane by a rotational isometry and so, without loss of generality, it is 
possible to restrict our attention to the study of equatorial geodesies only. 
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In what follows we denote the coordinate basis components by and the tangent vector 
to a curve by = dx"/dr. For timelike geodesies the parameter r can be made to coincide 
with the proper time and for null geodesies it only represents an affine parameter. Under 
the above circumstances, the geodesies take the following form (cf. IWaldl (119841 )): 



K = gapu'^u'' = -1 h'+ l + rV', (2) 



where 



f 1 for timelike geodesies, 
\ for null geodesies. 

In the derivation of the geodesic equation ([2]), there are two important constants of motion 
that must be taken into account. The first of them is 



where represents the static Killing vector and E' is a constant of motion. For timelike 
geodesies E represents the specific energy of a single particle following a given geodesic, 
relative to a static observer at infinity. 

The second constant of motion h is related to the rotational Killing field by the 
following relation: 



h:= g^ai^'^u'' = r^sm^9^. (4) 

dr 

Since we have chosen 6 = tt /2 without loss of generalisation, the previous equation takes the 
form 



h = r^^. (5) 
dr 

For timelike geodesies h is the specific angular momentum. The final equation for the 
geodesies is found by direct substitution of equations ([3]) and ([5]) into relation ([2]). From now 
on, we restrict the analysis to timelike geodesies only and so, the equation of motion takes 
the following form: 
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dry \ J \ ^ 

This equation shows that the radial motion of a geodesic is very similar to that of a 
unit mass particle of energy E"^ in ordinary one dimensional non-relativistic mechanics. The 
feature provided by general relativity in equation ([6]) is that, apart from a Newtonian grav- 
itational term —2M/r and the centrifugal barrier h? jr^^ there is a new attractive potential 
term —IMl^jr^^ that dominates over the centrifugal barrier for sufficiently small r. 

As it is done in the a i ialysis of the Keplerian orbit for Newtonian gravity (see for 



example iLandau &: Lifshitzl (jl994bl )). it is useful to consider r as a function of instead of 



r. Therefore, equation ([6]) takes the form 

dr \ ^ 2Mr3 



+ 2Mr + (E^ - l) . (7) 



Now, letting E"^ — 1 := 2£'tot, where i^tot is the total energy of the particle and u = r ^ 
equation ([7]) takes the final form 



duV _ o 9 2Mu 2Etot 



dip) ' 
Let us define u := Mv/h"^ so that the previous equation simplifies to 



dw^ ^ 



dcf> 
where 



av^ -v'^ + 2v + e, (9) 



mV 2E,.,h^ 



This equation determines the geometry of the geodesies on the invariant plane labelled 
by = 7r/2. In fact, this equation governs the geometry of the orbits described on the 
invariant plane due to the fact that the geometry of the geodesies is determined by the roots 
of the cubic equation 



f(v) = av^ -v'^ + 2v + e. 



(11) 
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The parameter a provides the difference between the general relativistic and the Newto- 
nian case. In fact, a — > in the Newtonian hmit. Finally, the eccentricity e of the Newtonian 
orbit is related to e through the relation 



1 + e. 



Accretion model 



The model first proposed by lUlrichI (119761 ) describes a non-relativistic steady accretion 
flow that considers a central object for which fluid particles fall onto it due to its gravita- 
tional potential. Their initial angular momentum at infinit y is considered small in such 



a way that this model is a small perturbation of iBondil (119521 ) 's spherical accretion model. 
The specific initial conditions far away from the origin combined with the assumption that 
radiative processes and viscosity play no important role on the flow, imply that the stream- 
lines have a parabolic shape. When fluid particles arrive at the equator they thermalise their 
velocity component normal to the equator. Since the angular momentum for a particular 
fluid particle is conserved, it follows that particles orbit about the central object once they 
reach the equator. The radius rdN o f the Newto r iian accretion disc, w here particles orbit 
about the central object, is given by (IUlrichlll976l : iMendoza et al.ll2004l ) 



rdN = hl,lM. 



(12) 



The velocity field and the density profiles are calculated by energy and mass conservation 
arguments. 



We consider now a general relativistic lUlrichI situation in which rotating fluid particles 
fall onto a central object that generates a Schwarzschild space-time. As described in sec- 
tion [H our analysis is well described by a ballistic approximation. The equation of motion 
for each fluid particle is thus described by relation ([9]). 

In order to get quantitative results it is important to establish the boundary conditions 
at infinity. The angular momentum is given by equation ([S]), so if a particle that falls onto 
the black hole has an initial velocity vq at an initial polar angle and the radial distance 
between the particle and the black hole is ro, then the angular momentum is given by 



2d(^ 



^01- = ^qIqVq sin 6*0, 
dr 



(13) 
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where 70 : = 



-1/2 



is the Lorentz factor for the velocity Vq. 



In addition, /iqo* is related to the angular momentum h perpendicular to the invariant plane 
through the relation 



h = /if>o* sin^n. 



(14) 



In the Newtonia n case, the specific angular momentum /iqo* converges to the value calculated 
bv luirichi Jl976h . 



With the above relations it is then possible to calculate the equation for a given fluid 
particle falling onto the central object. First of all, equation (ITT]) states that, if f{v) is a 
cubic polynomial in v, then either all of its roots are real or one of them is real and the two 
remaining are a complex-conjugate pair. The fact that the particle's energy is insufficient to 
permit its escape from the black hole's gravitational field requires that e < 0. This implies 
that the roots Vi, V2, and ^3 of f{v) are all real and satisfy the inequality fi < t'2 < ^3. 
Thus, /(f) can be written as 



f{v) =a{v- vi) {v2 - v) (f 3 - v) 



(15) 



Direct substitution of this relation in equation (Q yields the integration 



dw 



[{V2 - Vi) {Vs - fl)]^^^ J [{W^ - Wl) {W^ - Wi)]^^^ 



(16) 



where wf = 1/ (f2 — fi), W2 = 1/ (f3 — fi) and v = vi + w~'^. This elliptic in t egral can 
be cal culated in terms of Jacobi elliptic functions (see for example ICayleyl (jl96ll ): [Hancock 
( I1917I )) yielding the following result 



^ns \ w(v2- ViV^'^ 



{V3 - Vi)^^"^ 

The modulus k of the Jacobi elliptic function for this particular problem is given by 



;i7) 



V2 - Vi 



V3 - Vi 



With the aid of relation fll7p . the equation of the orbit is now obtained: 
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V = vi + {v2- vi) sn^l^ [a (^3 - ^^i)]^^^ }• (19) 

This is a general equation for the orbit. For the particular case we are interested in, it must 
resemble the orbit proposed by Ulrich when a = 0. Thus, the equation of the orbit must 
converge to a parabola in this limit. This is possible if and only if the eccentricity e = 1, 
which in turn implies e = 0. All these conditions mean that the roots of equation (fTOl) are 
given by 

l-(l-8a)'/' l + (l-8a)'/' , , 



and so, the modulus k of the Jacobi elliptic functions in equation (USl) takes the form 

,o l-(l-8a)i/2 



1 + (1 -8a)i/2- 

Note that the previous equations restrict the value of a in such a way that 



(21) 



< a < 1/8. (22) 

When a = 0, Ulrich solutions are obtained and the case a = 1/8 corresponds to the case for 
which the angular momentum h = 4M = 2rg reaches a minimum value. 

The orbit followed by a single particle falling onto a Schwarzschild black hole with the 
Ulrich prescription is then given by 



r 

where 

1/2 

n : = 

M M 



^ V2 snV/5, (23) 



^ (avsf^^ /l + (l-8a)'/M /^L 2 . 2. 



and p is the latus rectum of the generalised conic. Note that in the Newtonian limit, the 
length defined by equation flM|) converges to the radius of the Newtonian disc rdN as 
shown by relation f|T2|) . 
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Before using the equation of the orbit to find out the velocity field and the particle 



( Cavlev 


1961; 


Hancock 


1917) 



sn^^z, k) + cn^^z, k) = 1 
sn(z, k) — ^ sin(2;), cn(2;, k) cos(2;), dia{z, k) 1, as — > 0, 
d d 

— sn(2;, k) = cn{z, k) dn(2;, fc), —cn{z^ k) = —sn{z, k) dn(z, k). (25) 

OlZ (XZ 

The relativistic conic equation is obtained by direct substitution of these relations onto 
equation (J23i) . giving 

^ (26) 



V2 (1 - Cn2(y9/3) 



This orbit lies on the invariant plane 9 = 7r/2. We now obtain an equation of motion in 
terms of the polar coordinate 6 and the initial polar angle made by a particle when it 
starts falling onto the black hole. To do so, we note the fact that in order to recover the 
geometry of the spherical 3D space as « ^ it should be fulfilled thatQ 

2 cn^gp/j + cn^g/j-l 

= 2cn^go/3-l • ^''^ 

Since the invariant plane passes through the origin of coordinates, then the radial coor- 
dinate r remains the same if another plane is taken instead of the invariant one. Therefore, 
the angle Oq is the same as the one related to the value of the angular momentum of the 
particle at infinity (cf. equation f|T3|l ). Thus, the equation of the orbit is found by direct 
substitution of equations ( fT3l) and (l27l) into ( l26l) . and is given by 

^ n sin^go(2cn^go/3-l) 
In order to work with dimensionless variables, let us make the following transformations 



^ Ulrich ( 1976h showed that cosip — cos 9/ cos6'o using geometrical arguments. For the general relativistic 



limit, one is tempted to generalise this result to cn(p/3 = cn0f3/cn9of3. However, this very simple analogy 
does not reproduce the velocity and particle number density fields for the Newtonian limit. 
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r 



Vi_ 



{i = r, 9, v?), 



where 



da;* 



M 



no :-- 



Vk-=\ — 



n 
no 



n. 



1/2 



r* J 



In the previous relations, the mass accretion rate onto the black hole is represented by M. 
The velocity fk converges to the Keplerian velocity of a single particle orbiting about the 
central object in a cir cular orbit w hen a = 0. The particle number density no converges to 
the one calculated by iBondil (119521 ) in the Newtonian limit for the same null value of a. 



Under the above considerations, the equations for the streamlines r{6), the velocity field 
Vr, vo,Vy, and the proper particle number density n are given by 



sin^ ^0 (2cn2^o/5 - 1) 
V2 {cnWo(3 - cnW(3) ' 



o -1/2 n cnpe snpe dnf3e 1/2 , 
-2r ' (3 33^^ // [6, 00, V2, 13) 



V0 = r 



v^ = r 



n 



sin 9 

/2cn^6'o/5 — cn" 



sin^ 

sin 6^0 /f 2 (cn^6'o/3 - cn^6'/3) 



sin e V 2cn2^o/5 - 1 



1/2 



(29) 
(30) 
(31) 

(32) 
(33) 



where the functions /i and /2 are defined by the following relations: 



r Q 2 sin'e {2cn'eof3-l)-V2 sin^ (cn^gp/? - cn^g/?) 

°' ■ (2cn2^o/3 - 1) {{cnWol3 - cnW(3f + (2 (3cn(3e snpO dn(3ef}' 
f2 {e,eo,V2,f3) := (301(300 sn(3eo dn/^^p + {sin^p cos^p {2cn^eo(3 - l) - 
-2/?cn/?^p snpOo dnpOo sin^^p} /v2r. 

Equations fl29l) - fl33l) are the solutions to the problem of a rotating gas sphere onto a Schwarzschild 
black hol e, i.e. t hey r epresent a relativistic generalisation of the accretion model first pro- 
posed by lUlrichI (119761 ). 
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4. Convergence to known accretion models 



We have mentioned before (cf. section [3]) that the analytical solution must converge 
to the Ulrich accretion model when a ^ 0. In order to prove this, note that three very 
important conditions are fulfilled when a — > 0: (a) the modulus k of the Jacobi elliptic 
functions vanishes, (b) the root V2 —>■ 2, and (c) the parameter /3 — > 1/2. These conditions 
together with equation (125!) imply that relations (I29l)-(l33l) nat urallv converge to the non- 
relativistic Ulrich model (see for example iMendoza et al. that is: 



sin^ 9() 



1 — cos 6 cos 6o ' 
cost^ ^ 



1/2 1 + 



COS^r 



v^ = r 



cos — COS Q f cos Q \ 
sin Q \ cos / 

j^y2sin6'o A cos6'\^^^ 
sin Q \ cos Qq ) ' 



1/2 



p = r 



^3/2 1 + 



COS f 

COS 6q 



-1/2 



{l + 2r-iP2(cos^^o)} 



(34) 
(35) 
(36) 
(37) 
(38) 



where P2{x) is the Legendre second order polynomial given by P2{x) •= (Scos^x — 1) /2. 

On the other hand, if we consider a particular case for which the angular momentum of 
the fluid particles is null, then equations fl29|) - (l33|l converge to 



n 



2-l/2^-3/2_ 



(39) 



These equations describe a radial accretion r aodel onto a Schwarzschild black hole. They 
correspond to the model first constructed by iMichell (Il972l ) when pressure gradients in his 
calculations are negligible. 



5. The extreme hyperbolic model 

As mentioned in section [3l the parameter a reaches its maximum value when a = 1/8, 
which corresponds to a minimum angular momentum h = 2rg. In this limit the module k 
of the Jacobi elliptic functions is such that A; = 1, f2 = 4 and P = ^/S as can be seen from 
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equat ions fl20|) . flMl) and fl27j) . Also, when A; ^ 1, the following identities are valid (ILawden 
1989h : 



sn w — i> tanhw;, cnw — sechw, dnw— i>sechw. 
Using all these relations it follows that solutions fl29|) - fl33|) converge to 



(40) 



where 



/iH (d, 0o) 



Vg = r 



sin^ 00 ('2sech2^^o - 1 

r = — -p 

4(sech^^^n -sech2#^ 



_y,V2 sech^^g tanh^g 

-Y 

^,/2sech2^^o-sech2^^ , 



sin 9 



,/,sin^^o ( (sech^^ ^0 - sech^^^^ 



1/2 



n 



2sech2#^ 

^-3/2 gjj^ 



4 ^0 - 1 



2/^/' (9, 9o) /2H (^, 



(41) 

(42) 
(43) 

(44) 

(45) 



2H 



2 sin2 ^ (^2sech2^^o - l) - 4 sin^ (sech^^^o - sech^^^) 

v^sech^^^ tanh^^^)'!' 

^0-1 I - 



(2sech^^9o - l) \ {sech^^9o - sech^^^^)' + f ^: 



9,9q) := ^^sech^^^6'o tanh^^^o + s sin^o cos^^o I 2sech^^^i 
— ^sech^^p6'o tanh^p6'o sin^ 6*0 > /4r. 



This model does not formally represent a relativistic Ulrich solution, since the orbit 
followed by a particular fluid particle has a hyperbolic Newtonian counterpart. The solutions 
described by equations (l^-(H5ll ar e the e xact relativist ic solutions to the nu r nerica l problem 
discussed by iLee &: Ramirez-RuizI (l2006l ) who used a IPaczynsky fc Wiital (1l980l ) pseudo- 
Newtonian potential. 
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Discussion 



Ulrichl 's Newtonian accretion model predicts tlie existence of an accretion disc of radius 



TdN- This is a natural property of an accreting flow with rotation and has to be vahd in 
the relativistic case as well. In order to see the modifications that a full relativistic model 
imposes to the structure of the accretion disc, let us start by obse rving what happens to 
a fluid particle when it reaches the equator. First, in the lUlrichI accretion model, when 
any particle reaches the equator 9 = 7i/2, it does so at a radius r = h?/M according to 
the dimensional form of equation (!34|) . This corresponds to a stable circular orbit about the 
central object only in the case of a particle with azimuthal velocity that lies on the equatorial 
plane. For the relativistic model we have discussed s o far, if this were the case, then particles 
would arrive at the equator at a radius (IWaldlll984l ) 



r-circ = |^|l + (l-6a)^/^|, (46) 

which corresponds to the radius of stable circular orbits. However, when 9 = 7r/2, equa- 
tion (1291) implies that the value of r is very different from the one that would appear if a 
stable circular orbit is expected according to equation fH6|) . In fact, fluid particles arrive at 
a radius greater than rcirc- 

We can also discuss what happens to the radius of the disc for any a. This radius 
is obtained by taking a particle that arrives from a streamline just above the equator, i.e. 
^0 = 7r/2 — ?7, where the positive quantity ?7 -C 1. Figure [1] shows how this radius varies 
as a function of a. As it can be seen, the radius grows monotonically from the value 

when a = to infinity when a = 1/8. This behaviour strongly modifies the traditional 
view, particularly since the disc occupies all the equatorial plane in the extreme hyperbolic 
model. The fact that the disc radius diverges when a = 1/8 can be prooved directly using 
the results obtained in section [51 Indeed, evaluating equation fHTl) for 9 = 7r/2 and then 
taking the limit when 6*0 ^ 7r/2 it follows that r —>■ oo. 

The fact that the disc radius grows monotonically as a approaches the value 1/8 means 
that the density of the disc should be distributed in a more homogeneous way. Figure [2] 
shows density profiles evaluated in the equatorial plane ^ = 7r/2 as a function of the distance 
to the central object. In all cases the particle number density diverges in the origin because 
it represents a point of accumulated material. The case a = corresponds to the non- 
relativistic Ulrich model and apart from the divergence that the particle number density 
has at r = it also grows without limit at the radius of the disc raN- This is generally 
attributed to border effects that appear because the disc has been assumed to be thin (see 
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Fig. 1. — The figure shows a plot of the radius of the disc r' measured in arbitrary units, as 
a function of the parameter a. In the non-relativ istic ca s e, for which a = 0, the radius of the 
disc is exactly the same as the one predicted by lUlrichl (Il976l ). For the extreme hyperbolic 
model, when a — 1/8, the radius of the disc grows without limit. 
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e.g. lMendoza et al.ll2004j . and references therein). However, as Figure [2] shows, the divergence 
of the particle number density at the border of the disc disappears as soon as a moves away 
from a null value. Furthermore, it does so in such a way that the density of the disc varies 
very smoothly throughout the disc as a —>■ 1/8. 



The results of section [5] ca n be used to compare with th e pseudo-Newtonian iPaczynsky fc Wiita 



(ll980l ) approximation used by lLee fc Ramirez- Ruia (120061 ). Figure [3] shows a comparison be- 
tween the full relativistic solution with the pseudo-Newtonian approximation. It is clear 
from the images that the solution differs not only at small radii, near the Schwarzschild 
radius, but also at large scales. This is due to the strength of the gravitational field pro- 
duced by the central source, which makes particles approach the equator quite rapidly. For 
instance, near the event horizon there are fluid particles that appear to be swallowed by the 
hole when described by a pseudo-Newtonian potential. However, the complete relativistic 
solution shows that for this particular case some of those particles are not swallowed directly 
by the hole, but are injected to the accretion disc. 

The work presented in this article r epreseri t s a ge neral relativistic approach to the New- 
tonian accretion flow first proposed by lUlrichI (119761 ). The main features of the accretion 
flow are still valid with the important consequence that, the radius of the equatorial accre- 
tion disc grows from its Newtonian value for the Ulrich case up to infinity in the extreme 
hyperbolic situation, for which the angular momentum is twice the Schwarzschild radius. As 
a consequence, the particle iiumber density diverges on the border of the disc only for the 
Newtonian case described by lUlrichl . This is due to the fact that, when the radius of the disc 
grows, the particle number density on it rearranges in such a way that it smoothly softens 
as the extreme hyperbolic case is approached. Figures H] and [5] show streamlines and density 
isocontours for different values of the parameter a. 
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Fig. 2. — The plots represent different particle number densities n measured in units of no, 
as a function of the radial distance R (measured in units of r^) evaluated in the equator, i.e. 
for which the polar angle 9 = 7i/2. From bottom to top the models correspond to values 
a of 1/8, 10^^, 10^^, 10^'^, 10"^ and 10~^. All profiles diverge at the origin because of 
accumulated material at that point. The particle number density diverges in the Newtonian 
limit (for whic h q — » 0) at the bor der of the disc, which corresponds for that particular 
case to i? ^ 1 flMendoza et al.ll2004l ). However, this singularity disappears and softens the 
density profile in the disc as a — >• 1/8. 
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Fig. 3. — The figure shows a comparison between t he fully relativistic solution presented 
in this article (continuou s lines) with the New t onian iPaczynsky fc Wiital (119801 ) numerical 
approximations made by iLee fc Ramirez- RuizI (120061 ) (dotted lines). Distances in the plot 
are measured in units of the Schwarzschild radius. The plot is a projection at an azimuthal 
angle (p = const. The length R is the radial distance measured in the equator. In both 
cases, the streamlines were calculated in the extreme hyperbolic case for which a = 1/8, 
i.e. the specific angular momentum for a particular fluid particle is twice the Schwarzschild 
radius. Particles were considered to be uniformly rotating at a distance of 50 Schwarzschild 
radius measured from the origin. Both, the small and large scale panels show that the com- 
plete relativistic solu tion differs significantly from their calculations. The pseudo-Newtonian 



Paczynsky fc Wiital approximations were kindly provided by W. H. Lee. 
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Fig. 4. — Streamlines for values of the parameter a = 10"^, 0.12 from left to right are 
shown in the figure. Lengths are measured in units of the radius r^,. The equatorial radius 
is labelled by R. Th e case a = 10~^ is very close to the Newtonian one (see for example 



Mendoza et al.l (120041 ) ). This particular case shows that the streamlines are accumulated at 



R = 1, which corresponds to the Newtonian radius r^N- However, the right panel shows that 
as a approaches the value 1/8 the streamlines are not packed together any longer. 
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Fig. 5. — Particle number density isocontours for a = 10 ^, 0.05, 0.12 are shown in each 
diagram. The l eft pa nel roughly corresponds to the non-relativistic described by 



Mendoza et al.l (120041 ). All models show a density divergence at the origin. However, only 
the Newtonian case exhibits another divergence at the border of the disc R = 1. Lengths 
in the plot are measured in units of the radius r* and the density isocontours correspond to 
values of n/rio = 0.1, 0.6, 1.1, 1.6, 2.1, 2.6. 



